Curvature-induced anisotropic spin-orbit splitting in carbon nanotubes 
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We have theoretically explored the spin-orbit interaction in carbon nanotubes. We show that, 
besides the dependence on chirality and diameter, the effects of spin-orbit coupling are anisotropic: 
spin splitting is larger for the higher valence or the lower electron band depending on the specific 
tube. Different tube behaviors can be grouped in three families, according to the so-called chiral 
index. Curvature-induced changes in the orbital hybridization have a crucial role, and they are shown 
to be family-dependent. Our results explain recent experimental results which have evidenced the 
importance of spin-orbit effects in carbon nanotubes. 
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Improvements in the quality of carbon nanotubes 
(CNTs) have enabled the fabrication of quantum dots 
aiming at the realization of spintronics devices 
CNTs present a high Fermi velocity and a twofold orbital 
degeneracy originating from the topology of the honey- 
comb lattice. The unique fourfold degeneracy of CNTs 
energy states (spin plus orbital moment) has been ob- 
served in CNT quantum dots (QDs) by magnetic field 
spectroscopy measurements [HI and makes them particu- 
larly interesting since, besides the spin degree of freedom, 
they present the orbital moment to allow for quantum 
manipulation. In a recent experiment spin-orbit cou- 
pling has been directly observed in CNT as a splitting of 
the fourfold degeneracy of a single-electron energy level 
in ultra-clean QDs. This important finding seems to be 
in contradiction with the interpretation of earlier experi- 
ments in defect-free CNTs, from which independent spin 
and orbital symmetries and electron-hole symmetry have 
been deduced Besides showing the importance of 
spin-orbit effects in carbon nanotubes, Kuemmeth et al. 
[6| point out an unexplained anisotropic splitting of elec- 
tron and holes in carbon nanotube quantum dots, which 
deserves further exploration. 

On theoretical grounds, spin-orbit interaction (SOI) 
has been investigated on CNTs by deriving an effec- 
tive mass Hamiltonian including a weak SOI in carbon 
orbitals to the lowest order in perturbation theory Q. 
Band splitting was found considering surface curvature 
effects [91, as well as in the electron spin resonance spec- 
tra of achiral CNTs derived by low-energy theory [lO|. 
In an earlier work, we showed that the inclusion of the 
full lattice symmetry is essential for deriving spin-orbit 
(SO) effects in CNTs Employing an empirical tight- 
binding model, we demonstrated an intrinsic symmetry 
dependence of SOI effects. As confirmed by recent ex- 
perimental results Q, we showed that, in the absence 
of a magnetic field, CNTs present spin-orbit split bands 
at the Fermi level. In addition, SOI induces zero-field 
spin splitting in chiral CNTs, while Kramers theorem on 
time-reversal symmetry alongside the inversion symme- 



try preserve the spin-degeneracy in achiral — i.e., {n,n) 
armchair and (n, 0) zigzag — nanotubes [13] ■ More re- 
cent works [13I I have indicated the importance of cur- 
vature in the SOI effects investigated with a continuum 
model by perturbation theory, reporting the appearance 
of a gap and spin-splitting in the CNT band spectrum 
due to SOI. Here we show that, besides the dependence 
on the diameter and chirality, SOI effects in CNTs ex- 
hibit an electron-hole anisotropy which is specific to the 
tube. Different nanotube behaviors with respect to SOI 
are grouped into three families, which have also arisen 
with respect to other electronic properties [14|. 

We have performed electronic structure calculations 
using an empirical tight-binding Hamiltonian including 
a four orbital sp^ basis set 15|, ll6| . The CNT unit cell is 
formed by rolling up a portion of a graphene sheet; thus, 
the actual discrete nature of the lattice and curvature 
effects are taken into account. The atomic spin-orbit in- 
teraction term included in the Hamiltonian is given by 
Hso — -^L • S, where A is a renormalized atomic SO cou- 
pling constant and L and S stand for the orbital and spin 
angular momentum of the electron, respectively. The 
spin quantization direction has been chosen parallel to 
the carbon nanotube axis. Different estimations of A 
have been done [1, [3], always assuming a very small 
value, considerably reduced in graphite/graphene and 
CNTs with respect to the value for atomic C (« 12 meV). 
However, the recent experiments of Ref. ^] point out an 
enhancement of the role of SOI in CNTs with respect to 
that of graphene; thus, the exact value of the SO cou- 
pling parameter is still under discussion. In the present 
work, SO-induced energy splittings are given relative to 
the strength of the SOI, and only in band structure re- 
sults we have chosen an artificially large value of A, for 
the sake of clarity in the Figures. 

We focus first on achiral zigzag nanotubes with a chi- 
ral angle = 0°, where curvature effects are expected 
to be largest 17, 3]. Zigzag (n,0) tubes can be classi- 
fied as primary metals 18| if n = 3q, q being an integer, 
or semiconducting, if n = 3q ± 1. Without curvature 
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FIG. 1: Band structures calculated around the Fermi level 
for three zigzag CNTs: (8,0) tube from the (3g — 1, 0) family, 
(9,0) primary metallic nanotube from the (3g, 0) family, and 
(10,0) tube from the (3g + 1,0) series. Band structures have 
been calculated including the SOI term with A — 0.2 eV. 



effects, primary metallic zigzag tubes have a band cross- 
ing at the Fermi level; this occurs at the center of the 
Brillouin zone (BZ). The bands crossing at F are four- 
fold degenerate, and inclusion of curvature opens a small 
gap. Zigzag tubes with n — ± 1 are semiconduct- 
ing with the bandgap at F; the top valence and bottom 
conduction bands are also fourfold degenerate. Inclusion 
of SO interaction partially removes the band degeneracy 
at Ef, although the split bands remain spin-degenerate. 
The energy splittings induced by the SO interaction term 
are different for the highest valence band (VB) and the 
lowest conduction band (CB): in two of the families, for 
CNTs (3q, 0) and (3g — 1, 0), the splitting is larger for the 
VB, whereas for NTs {3q + 1, 0) the splitting is larger for 
the CB. In Figure [T] we show one particular example of 
each of the three zigzag families, namely the (8,0), (9,0) 
and (10,0) CNTs. Although our analysis concentrates 
on those bands closest to the Fermi level, it can be seen 
that the splitting of all other bands is also anisotropic 
and specific to each band. In particular, note that the 
second CB and VB show opposite behaviors to those of 
the bottom CB and top VB bands (see below). 

In Figure [2] the SO energy splitting of the top valence 
and bottom conduction bands is represented as a func- 
tion of the tube diameter for the three zigzag CNT se- 
ries, showing the family behavior described before. The 
curvature effect is clearly seen: the SO energy splitting 
is much higher for CNTs with the smallest diameters; 
it can be noted as well how the differences among the 
three series decreases with increasing diameter. We have 
checked that the splitting is not negligible for larger NTs: 
for example, the (80,0) tube, with a diameter of 62.64 A, 
roughly the same size of the NT measured in '6^|, has a 
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FIG. 2: (Color online) Normalized spin-orbit splitting of the 
top valence and bottom conduction bands as a function of the 
tube diameter for the three zigzag families. 



top VB normalized splitting of 0.04. An analogous fam- 
ily behavior of zigzag CNTs has been observed in other 
physical properties, such as the band gaps H, 2^, 21 1, 
and it has been related to the trigonal warping effect. 
The three families can be understood by resorting to the 
zone folding approach. Within this approximation, the 
energy bands of a general (n, m) CNT are given by impos- 
ing periodic boundary conditions to the graphene sheet: 
if the lines of allowed k vectors touch the Fermi points 
{K or K') of the hexagonal graphene BZ, the nanotube 
is a metal. This occurs when n — m = 3g, with q integer; 
otherwise it is a semiconductor. The two semiconduc- 
tor families, n — m = ± 1, correspond to nanotubes 
for which the quantization lines yielding the energy gap 
are at opposite sides of the Fermi point K (or K'). Thus, 
CNTs can be classified into three families: n~m = iq+v^ 
with V = Q, ±1 being the so-called chiral index, hereafter 
referred as family index [2^]. Denoting as K the BZ spe- 
cial point with coordinates (|f,0), the v — ~1 has the 
closest quantization line yielding the gap to its left, over 
the TK line, whereas the v = +1 family has it to the 
right, over the KM line. Of course, if we choose the K' 
Fermi point to classify the nanotube families, the relative 
positions (right/left) of the quantization lines yielding 
the gap are reversed with respect to K . For simplicity, 
henceforth we will refer our discussion to the K point 
with coordinates given above. 

In Fig. [3] we illustrate this general classification for the 
case of zigzag CNTs, with two examples belonging to dif- 
ferent series. The allowed quantization lines of the (7,0) 
and (8,0) tubes closer to K are shown over the graphene 
TT-band structure energy contour plot. The trigonal warp- 
ing effect is clearly seen in these contours: quantization 
lines at opposite sides of K correspond to graphene bands 
with appreciably different slopes, although they are at 
similar distances. As the tt graphene bands crossing at 
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FIG. 3: (Color online) Contour plot of the graphene band- 
structure calculated around the Fermi K point. The quan- 
tization lines corresponding to the (7,0) (green, dashed) and 
(8,0) (blue, dotted) CNTs closer to this point are also shown. 



K have different symmetry, the character of the conduc- 
tion band changes from antibonding along the TK line 
to bonding along the KM line. This explains why tubes 
belonging to = ±1 families show an opposite behavior 
in their conduction (or valence) bands: the two bands 
closest to the Fermi energy have a different symmetry. 
In addition, it clarifies why in a given tube, as described 
above, the second VB (CB) has a reverse behavior with 
respect to the first VB (bottom CB): these second bands 
arise from quantization lines at opposite sides of K with 
respect to the first. Furthermore, the similar behavior of 
the metallic and the semiconductor v = —1 families can 
be explained by noticing that curvature effects shift the 
Fermi K point in Fig. [3] to the right, in such a way that 
the allowed quantization lines of the metallic tubes closer 
to K fall onto the same side as those corresponding to 
the 1/ = -1 family [I^. 

On the other hand, curvature changes the hybridiza- 
tion of the orbitals. The linear graphene bands crossing 
at K are of a pure tt character around Ep] however, re- 
hybridization of a and tt orbitals is very important for 
small-radii nanotubes, and non-negligible in general. In 
order to show the a-ir hybridization of the zigzag nan- 
otubes we have calculated the contribution of each of the 
four orbitals forming the basis set for the valence and 
conduction band states closest in energy to the Fermi 
level at F; the results summed over all the atoms of 
the unit cell are shown in Table HI NTs belonging to 
the v = +1 family have a larger contribution of the a 
orbitals — that is, a larger ct-tt rehybridization — in the 
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conduction band than in the valence band; however, in 
tubes of the v = 0,-1 families the a orbital density is 
larger in the valence band. Spin-orbit effects are more 
important for bands with larger curvature-induced rehy- 
bridization: therefore, curvature effects are responsible 
for the observed electron-hole anisotropic SO splitting. 
The different mixing of tt and a orbitals, due to the cur- 
vature of the tubes, was shown to affect the energy gaps 
of semiconducting nanotubes [l3| with a similar family 
dependence. Here we have shown that it also influences 
the value of spin-orbit splitting in CNTs. 



TABLE I: Electronic densities 
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We consider now chiral NTs (n,™), with 7i ^ m ^ 
that do not have an inversion center. As mentioned 
above, curvature effects induce a shift of the Fermi wave 
vector kp, opening a small gap at the Fermi energy in the 
primary metallic chiral CNTs, (m — n) = 3q. But, in con- 
trast to the results shown for achiral zigzag NTs, in both 
metal and semiconductor chiral NTs the SO interaction 
lifts all degeneracies. 

The calculated energy splittings for these tubes follow 
the same behavior as those obtained for achiral zigzag 
CNTs: for tubes with v = {n — m ^ 3q) and — —1 
(n—m = 3(?— 1) the energy splitting is larger for the high- 
est VB, while for tubes with i/ = +1 {n — m = 3q + 1) 
the splitting is larger for the lowest CB. As an example 
for chiral tubes, the band structures calculated includ- 
ing the SOI for three particular tubes, (6,4), (9,3), and 
(8,4) belonging to each of the three families, are shown 
in Fig. [4] around the Fermi level. However, chirality has 
an important effect in SOI: this is illustrated in Fig. [5l 
where the band splittings vs. diameter for chiral tubes 
belonging to the three families are depicted. The symbol 
color indicates the NT chiral angle. As a guide to the 
eye, the zigzag metallic i/ = results are also shown. It 
can be seen how chiral tubes follow the family behavior 
already described; notwithstanding, chirality effects in- 
troduce deviations from the zigzag monotonic behavior, 
due to the different orientation of the relevant quantiza- 
tion lines with respect to the K point. In general, the 
higher the chiral angle, the larger the deviation from the 
behavior of zigzag nanotubes. 

In summary, we have shown that spin orbit effects in 
carbon nanotubes is anisotropic: The energy splittings 
induced by the SO interaction term are different for the 
highest valence band and the lowest conduction band de- 
pending on the tube family; the magnitude of the SO 
splitting correlates with the a-n hybridization induced 
by curvature. These dissimilar a-n hybridizations of the 
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FIG. 4: Band structures calculated with SO interaction 
around the Fermi level for three chiral CNTs: the (6,4) tube 
from the = —1 family, the (9,3) primary metallic nanotube, 
and the (8,4) tube from the v = +1 family. 
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FIG. 5: (Color) Normalized energy band splittings for chiral 
CNTs versus diameter. Squares stand for i/ = +1 family, 
triangles for u ~ 1 and circles for = ; full (open) indicate 
the top VB (bottom CB) splittings. Symbol color indicates 
the nanotube chiral angle. Full (dashed) line: Bottom CB 
(top VB) splittings for metallic zigzag tubes. 



valence and conduction band states, which in turn de- 



pend on the position of the quantization hnes yielding 
the gap with respect to the graphene K point, are the 
reason for the experimentally observed anisotropy in SO 
splitting in absence of external fields. 
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